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i 16. Abstract 

i 


Using a 3-compancnt sinusoidal tit of ^ versus sunspot cycle number (where Is the smoothed sunspot number nt 

eyclo maximum) for cycles 8 through 20, onsldered to bo the most reliably known cycles, values of R frIAX are projected for cycles 

21 and 22. For cycle 21, the fit project, an value of 160,5, which compares very well to the actually observed value of 

164.6, For cycle 22, it projects an value near 110, Results of linear regression analysis, based on comparisons to observed 

valuos of arc presented which allow for a determination of a number of cycle-related parameters, Including ASC, DES, MIN- 

MIN PERIOD, st, O pE ASC* 011,1 SI,ot ’ E DES (defined in the text). It is noted that sunspot cycles can be divided, according 

to their mlnlmum-to-mlntmum cycle duration, into two distinct groupings: (1) short-period cycies having a mean cycle duration of 

about 122 months (ranging between 116 and' 126 months) and (2) long-period cycles having# mean cycle duration of about 140 months 
(ranging between 134 and 149 months). Short-poriod cycles often arc observed as high-valued Ii, i!A j, cycles <4 of 6, with all cycles 

huvingH MAX > 142 being short-period cycles, 3 of 3), while long-period cycles often are observed as low-valued cycles 

(5 of 7). This, paradigm relating cycle amplitude to period (obsorvablo in 0 of the 13 modern cycles) suggests that cycle 21, being 
n high-valued cycle (and with "^j AX > 142), may also be a cycle of short duration. If so, then cycle 21 will end about 

August 1936. In contrast to this, however, a plot of MIN-MIN PERIOD versus sunspot cycle numbor, based on cycles 8 through 20, 
suggests that cycle 21 may be a cycle of long duration. Consequently, occurrence for cycle 22 would be delayed until about 

Fobruary 1938, Since cycle 22 Is projected to bo a low-valued K M ax cycle with a value very near the observed value for cycle 20, 

tho cycle Is expected to bo of long duration. Consequently, cycle 22 should peak about October 1990 and end about April 1998, If 
cycle 21 Is a short-period cycle ending about August 1906, or cycle 22 should peak about April 1992 and end about October 1993, 

if cycle 21 la a long-period cycle. [Two additional, but less likely, predictions must also bo noted. Those include: (1) if cycle 21 

1b a Ipng-poriod cycle and cycle 22 a short-period cycle, then cycle 22 should peak about October 1991 and end about February 
1098; and (2) If cydu 21 Is a short-period cycle and cycle 22 e short-period cycle, then cycle 22 should peak about April 1990 and 
end about August 1996. J Using SLOPE A g^ and 6LOPE DE g (based on projected and measured SLOPE Ag( , valuas) and a "transient" 

curve-fit scheme, smoothed sunspot numbers arc projected for cycle 21 and compared to observed values. This approach is adequate 
to describe tho gross behavior of smoothed sunspot number for eyclo 21. Following tho same approach, values of smoothed sunspot 
number arc projected for cycle 22. An abbreviated version [21J of this report is to be published in the Proceedings of the Solar 
Terrestrial Prediction Workshop (13-22 Juno 1984, Moudon, France, P. A. Simon, Workshop Coordinator), 
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TECHNICAL MEMORANDUM 
ESTIMATING SUNSPOT NUMBER 
I. INTRODUCTION 


The cyclic nature of sunspots has been a recognized solar phenomenon since 
about 1850, following the introduction by Rudolf Wolf in 1848 of the now well-known 
"relative sunspot number R" and the announcement by Heinrich Schwabe in 1843 of 
an apparent 10-year periodicity in sunspot observations made between 1826 and 1843. 
Until very recently, sunspot number had been routinely and systematically measured, 
following Wolf’s numerical method, by the Swiss Federal Observatory in Zurich, 
Switzerland and its two branch stations in Arosa and Locarno; it was denoted R^. 

Beginning in January 1981, an "international" sunspot number, denoted Rj, has been 

measured by the Sunspot Index Data Center in Brussels, Belgium and replaced R^. 

The highest quality sunspot number data are regarded to be those which have been 
reported since about 1850 (Eddy [1]). Consequently, sunspot cycles dating back to 
Schwabe's original sunspot cycle, denoted sunspot cycle number (SCN) 8, are con- 
sidered to be the most reliable and, thus, form the basis for a determination of 
statistical properties of sunspot cycles and for estimating the variation of Bunspot 
number with time. The last completed cycle is SCN 20 which began in October 1964, 
peaked in November 1968, and ended in June 1976. SCN 21 began in June 1976, 
peaked in December 1979, and is now in decline. SCN 22, the next cycle, is expected 
to begin in the late 1980's (either late 1986 to early 1987 or late 1987 to early 1988, 
dependent upon cycle 21 being a short-period or long-period cycle, respectively; see 
Section III.D). 


On the basis of smoothed sunspot number, denoted R^ 3 (also called "monthly 

moving average," "13-month running mean," etc.), for cycles 8 through 20, the 
following are examined; (1) sunspot cycle statistics, (2) an empirical curve-fit for 
^MAX versus ®CN where Rjyj^ x * s maximum R^ 3 value for a cycle marking cycle 


maximum , 


(3) linear regression equations based on 


H max . (4) the inferred bi-modal 


distribution of sunspot cycles based on cycle duration, (5) a transient curve-fit 
scheme, and (6) application of these techniques to cycles 21 and 22. Section II 
describes the data base upon which this study is founded. Section III discusses the 
above topics in the stated order, and Section IV gives the conclusions. Much of this 
work has been extracted from previous studies which have been published elsewhere 
(e.g., Teuber, et al. [2,3] and Wilson [4,5,6], 


II. APPROACH 


The Wolf relative sunspot number R is defined as 


R = k(10g + f) , 


( 1 ) 


where f is the total number of sunspots observed regardless of size, g is the number 
of observed sunspot groups, and k is a normalization parameter which varies from 
observatory to observatory to bring counts into agreement by accounting for telescope 
size, atmosphex’ic opacity, etc. Sunspot counts are made daily and then averaged at 
the end of the month to obtain a monthly mean sunspot number. As stated in the 
introduction, monthly mean sunspot numbers of historical importance are those made 
and collated by the Swiss Federal Observatory located in Zurich, Switzerland, denoted 
R z . Waldmeier [7] has compiled values' of R z for the period 1610 to 1960. Values of 

R z (now Rj) for the years since 1960 (actually, since 1944) are published monthly in 

Solar Geophysical Data— Prompt Reports (NO A A /Environment Research Laboratory, 
Eoulder, Colorado USA). 

For statistical comparisons, a smoothed sunspot number R 13 (also denoted in 
Waldmeier [7]) has come into use. It is defined as 


R +6 + R -6 + 2 Z R i 


R 13 = 


+5 

z 

i=-5 


24 


( 2 ) 


where R + g is the monthly mean sunspot number 6 months ahead of the month of 
interest, R_g is the monthly mean sunspot number 6 months behind the month of 
+5 

interest, and ^ R. is the sum of the monthly mean sunspot numbers 5 months 
i=-5 

either side and including the month of interest. In this study, use is_made of 
smoothed sunspot number. Wilson [5] has compiled values of R„ and R.n for cycles 
8 through 20 and part of cycle 21. L 


III. RESULTS AND DISCUSSION 


A. Sunspot Cycle Statistics 

Table 1 summarizes specific cycle-related parameters for cycles 8 through 20. 

It also gives parametric mean values and standard deviations for this group of cycles, 
and values for the parameters for cycle 21 when they are known. On the left, each 
cycle is identified by its sunspot cycle number (SCN). Adjacent to this and moving 
to the right are the parameters of interest; (1) the cycle minimum occurrence date 
(in month and year),_(2) the cycle maximum occurrence date (in month and year), 

(3) R MI jg (i.e., the R 13 value at cycle minimum occurrence), (4) E MAX (i.e., the 

R^g value at cycle maximum occurrence), (5) ASC (i.e., the ascent period or time 

in months from cycle minimum occurrence to cycle maximum occurrence), (6) DES (i.e., 
the descent period or time in months from cycle maximum occurrence to subsequent 
cycle minimum occurrence), (7) MIN-MIN PERIOD (i.e., the cycle duration or time in 
months from cycle minimum occurrence to subsequent cycle minimum occurrence, 
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TABLE 1. SUMMARY OF 3UNSPOT CYCLE INFORMATION FOR 

THE MODERN CYCLES 


SUNSPOT CYCLE CYCLE 

CYCLE 





MIN-MIN 




NUMBER (SCN) 

MINIMUM 

MAXIMUM 

»MIN "MAX 

ASC 

DES 

PERIOD 

_"MEAN SLOPEasC 

slope des 

8 

NOV 

1833 

MAR 1837 

7.3 

146.9 

40 

76 

116 

66.9 

3.490 

-1.795 

9 

JUL 

1843 

FEB 1848 

10.5 

132.0 

65 

94 

149 

56.9 

2,209 

-1.370 

10 

DEC 

1856 

FEB 1860 

3.2 

97.9 

60 

85 

135 

48.2 

1,894 

-1.091 

11 

MAR 

1867 

AUG 1870 

6.2 

140.5 

41 

100 

141 

52.9 

3,300 

-1.383 

12 

DEC 

1878 

DEC 1883 

2.2 

74.6 

60 

74 

134 

33.8 

1.207 

-0.941 

13 

FEB 

1890 

JAN 1894 

5.0 

87.9 

47 

96 

143 

38.5 

1.764 

-0.888 

14 

JAN 

1902 

FEB 1906 

2.7 

64.2 

49 

89 

138 

32.1 

1.256 

-0.704 

15 

JUL 

1913 

AUG 1917 

1.6 

105.4 

49 

71 

120 

44.0 

2.120 

-1.406 

16 

JUL 

1923 

APR 1928 

6.6 

78.1 

57 

65 

122 

40.3 

1.272 

-1.148 

17 

SEP 

1933 

APR 1937 

3.5 

119.2 

43 

82 

125 

57.2 

2.691 

-1.360 

18 

FEB 

1944 

MAY 1947 

7.7 

151.8 

39 

83 

122 

74.0 

3.695 

-1.788 

19 

APR 

1954 

MAR 1958 

3.4 

201.3 

47 

79 

126 

90.5 

4.211 

-2.427 

20 

OCT 

1964 

NOV 1968 

9.6 

110.6 

49 

91 

140 

59.9 

2.061 

-1.081 

21 

JUN 

1976 

DEC 1979 

12.2 

164.5 

42 

- 

— 

- 

3.626 

— 

MEAN VALUES (CYCLES 8-20) 

5.2 

116.2 

48.1 

83.5 

131.6 

53.4 

2.398 

-1.337 

STANDARD DEVIATIONS (CYCLES 8-20) 

2.7 

36.7 

6.2 

10.0 

10.0 

16,2 

0.958 

0.441 


numerically equal to ASC + DES), (8) ^ MEAN (i.e., the mean R 13 value calculated 

over the entire cycle)* (9) SLOPE A g C [i.e., the line slope on the ascent side from 

Kjvjjn ^ * computed as SLOPE^g^ — [Rj^ A jj.(SCN) *‘^j^j^(SCN)] /ASC(SCN)) , and 

(10) SLOPE DEg [i.e., the line slope on the descent side from R MAX to subsequent 

cycle E min , computed as SLOPE DEg = [K MIN (SCN + 1) - E MAX (SCN)] /DES(SCN)] . 

Clearly, sunspot cycles, on average, range in smoothed sunspot number from aboiit 
5.2 ± 2.7 (i.e., mean value ± 1 standard deviation unit) at to about 116.2 ± 

36.7 at ^ MAX and have a cycle duration of about 132 + 10 months, with the descent 

period being about 1.8 times longer than the ascent period. A superposed-epoch 
analysis, based on mean R^ 3 values computed as a function of time t from 

occurrence for cycles 8 through 20, yields parametric values close to those summarized 
in Table 1. One minor difference is that ^ MAX is slightly reduced to about 106.9 ± 

36.0. In terms of range, from the lowest observed parametric value to the highest 
observed value for cycles 8 through 20, it is seen that has varied from 1.5 to 

10.5, E max from 64.2 to 201.3, ASC from 39 to 60, DES from 65 to 100, and MIN-MIN 

PERIOD from 116 to 149. 


B. R max Versus SCN, An Empirical Curve-Fit 

In Figure 1, observed R MAX values are plotted against their respective SCN. 

For cycles 8 through 14, a downward trend is suggested, while an upward trend is 
suggested for cycles 14 through 20. For the whole set — i.e. , cycles 8 through 20 — 
a time series regression reveals an upward trend, with R„ AX being correlated posi- 
tively with SCN approximately as IAA 
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(3) 


**MAX " 90 + 2 SCN 
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The Pearson correlation coefficient r approximately equals 0.2 and the standard error 
of estimate S yx , corrected for small sample size, is about 40. 


| Rmax vs. SCN 
Is RmAX" 9 ° + 2SCN 



Figure 1. Linear regression fit of the variation of Kjy, . ^ 
versus SCN for modern sunspot cycles. 

Eddy [1,8-10] has done considerable research on the historical record of sun- 
spot data and solar activity. His studies, based primarily on naked-eye and early 
telescopic sunspot observations, fossil radiocarbon levels in tree rings, and auroral 
reports, have suggested possibly strong climatic associations in the Sun-Earth record; 
in particular, he has found evidence for the occurrence of prolonged "sunspot minima" 
(e.g., the Sp6rer Minimum — A.D, 1400-1510, and the Maunder Minimum — A.D. 
1645-1715) and a prolonged "sunspot maximum" (the so-called Medieval Maximum or 
Grand Maximum — ca., A.D. 1120-1280). If it is assumed that such sunspot minima 
and maximum have indeed happened and that they reflect periods when sunspot 
number was relatively lower and higher, respectively, equation (3) can be approxi- 
mated by equation (4), given below as 


^MAX = 12 ® + 25 s * n 


(4) 


Equation (4) is of the form y = a Q + a.^ sin 0-^Cx), where y is as before, and 

0^(x) is the phase of an SCN, within a 90-cycle periodicity. The coefficient aQ was 

selected have a value approximately equal to the value given by equation (3) for 
SCN 14, the midpoint of the modern data set. The coefficient a.^ represents a "best- 

fit" (to the sum of the squares of the residuals) for the given Sq and values of 0j(x) 
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contained in Table 2 , The 90-cycle periodicity of w ^ cJl * s used to approxi- 

mate the linear regression equation ^equation (3)], iB based on Eddy's suggestion 
that a Grand Maximum really existed, having a peak value at SON ^ -53, and that 
the Maunder Minimum likewise existed, The time difference between Grand Maximum 
and Maunder Minimum is about 500 years or 45 sunspot cycles, of 11 years ave rage 
duration; hence, the period of 90 cycles results. At SON :k -53, 0-^Cx) = 90 degrees 

is defined so that * s a * mQX ^ mum value. The 90-cycle periodicity causes each 

cycle to increment by 4 degrees; thus, for cycles 8 through 20, G-jfx) varies between 
334 and 22 degrees, respectively, (See Wilson [6] for additional comments.) 

TABLE 2. PHASE INFORMATION FOR THE 3-COMPONENT 

SINUSOIDAL FIT 


3 

J 

K 

? 

K 

SCN 

s 

5 

* 

0, IX) 

3 

$ 

0 2 IX) 

HI 

« 

0 3 IX) 

Umax <«>mp)1 

Umax < obs > 

a 

334 

90.0 

270 

129.0 

146.9 

9 

338 

122.7 

90 

155.1 

132.0 

to 

342 

165.5 

270 

118.8 

97.9 

tl 

346 

188.2 

90 

124.0 

140.5 

12 

350 

220.9 

270 

77.7 

74.6 

13 

35A 

253.6 

90 

98.8 

87.9 

14 

358 

286.4 

270 

70,6 

64.2 

IS 

2 

319,1 

90 

113.0 

105.4 

16 

6 

351,6 

270 

102.6 

78.1 

17 

10 

24.5 

90 

153,9 

119.2 

18 

14 

57.3 

270 

140.5 

151.8 

19 

18 

90.0 

90 

177,7 

201.3 

20 

22 

122.7 

270 

143,8 

110.6 

21 

26 

155.5 

90 

(160.5) 

164.5 

22 

30 

188.2 

270 

012.5) 

- 


1 *WxtC°**t » 1» * n tk I ,(K) ♦ X Mfldj (XI ♦ IS tin fljOO Ito. ISII 


VARIANCE- - M.Z * I* 

•*• t»» 

VARIANCE OX MIAN • ISM,* * 10.6 REDUCES VARIANCE »V 70* 


Returning to Figure 1, it is reiterated that there is an apparent downward trend 
in K MAX values for SCN 8 through 14 and an upward trend for SCN 14 through 20. 

Closer inspection reveals that, in addition to this variation, there is an "up -down- up" 
signature, especially between SCN 10 and 17. Together, these two observations 
suggest that the versus SCN plot might be better fitted using a 3-component 

sinusoidal curve, one component being a 90-cycle periodicity, mentioned above, a 
second component being an 11-cycle periodicity , and a third component being a 2-cycle 
periodicity. Such an analysis yields the relation 


5 


(5) 


RjyjAX “ ^0 + 33 s * n 0 1 ^ X ^ + 33 B * n + * 5 s * n 0 3 ^ x ) * 


where values for (^(x), e^Cx), and 6g(x) ai’e contained in Table 2. It is seen, then, 

that for any modern SCN, an approximate value for K MAX can be deduced, which is 

based on a fit that has a smaller variance than that for the mean. As compared to 
the mean, the total variance using equation (5) has been reduced by about 70 per- 
cent; thereby, achieving a standard deviation s of about 20. 


Figure 2 compares computed and observed values of K MAX for cycles 8 through 

21, where the computed R^x values are calculated using equation (5). The two 

curves behave similarly and all of the observed & MA x values He within 1.5 s units 

(i.e., ±30) of the computed R MAX values. It should be emphasized that the fit is 

based strictly on a fit of SCN 8 through 20, and, as such, it may be of little or no 
predictive value for future cycles. Its worth becomes apparent only after it has 
successfully predicted several successive cycles. It is noted, however, that equation 
(5) predicts SCN 21 to have an ^ MAX value of about ICO, which is remarkably close 


to its observed ^ MAX value, now known to have been about 105, a prediction accurate 


to within a few percent. 



Figure 2. Empirical fit using a 3-component sinusoidal 
function (R MAX versus SCN). 


C. Linear Regression Equations Based on R MAx 


Now that a reasonable estimate for Hjyj AX is possible, from equation (5), other 
cycle-related parameters can be estimated using results of linear regression analysis, 
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based on observed values of as the independent variable (Wilson [5,6]). 

Figures 3 and 4 depict regression equations for the parameters listed in Table 1, 
based on a comparison to Ri&u re 3 illustrates the variation of SLOPE A gQ, 


MEAN’ 


and with it MAX . Both SLOPE Agc and R^g^N show strong "positive" 


correlation, having a Pearson correlation coefficient r greater than 0.9. On the other 
hand, although correlated in a positive sense with K MAX , has an r of only 


about 0.3, indicating that the correlation is, at best, weak. Figure 4 plots MIN-MIN 
PERIOD, DES, and ASC against H^ x , ASC shows the strongest correlation of these 

three parameters, being negatively correlated with and having an r of about 

-0.(5. The rather strong correlations of SLOPE ASC and ASC with R MAX > while not 

independent, amply confirm that cycles* with high maximum sunspot numbers get their 
maxima more swiftly, rather than by sustaining their growth rates longer than cycles 
with lower maxima. Figure 5 plots SLOPEj^g against SLOPE A g^ and shows the strong 

negative correlation (r ^ -0.9) expected from the earlier, individual regressions 
against R MAX * 





Figure 3. Linear regression fits of 
SLOPE asc , n MEAN , and R MIN 

against & MAX . 


Figure 4. Linear regression fits of 
MIN-MIN PERIOD, DES, and 
ASC against E MAx . 



Figure 5. Linear regression fit of SLOPE DES against SLOPE ASC , 
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D. Inferred Bi-Modal Distribution of Sunspot Cyc9Js^^ QUALITY 
Based on Cycle Duration 

From Table 1, it is discerned that, although sunspot cycles average about 132 
months in duration, there have been no sunspot cycles in the entire set of SCN 8 
through 20 which have a cycle duration in the range 127 to 133 months. All cycles 
have either a minimum-to-minimum period shorter than 127 months or longer than 133 
months (Wilson [5] and Wilson et al. [11]). This bifurcation of cycle duration is 
readily seen in Table 3. Cycles 8 and 15 through 19 have cycle duration between 116 
and 126 months and cycles 9 through 14 and 20 have cycle duration between 134 and 
149 months. Excluding cycle 9 which has the longest cycle duration on record, being 
149 months, the range is between 134 and 143 months. Thus, sunspot cycles appear 
to be distributed by cycle duration into two distinct groupings: (1) a short-period 

cycle group, having a mean period of about 122 months with a standard deviation of 
about 3 months, and (2) a long-period cycle group, having a mean equal to about 
140 months and a standard deviation equal to about 5 months; excluding SCN 9 
reduces the mean for the long-period cycle group to 138.5 with a standard deviation 
of 3.2. A summary is provided in Table 4. 


TABLE 3. DISTRIBUTION OF CYCLE DURATION 

CYCLE DURATION FREQUENCY OF 
INTERVAL (MONTHS! OCCURRENCE SCN'S 


116 - 121 
122 - 127 
128 - 133 
134 - 139 
140 - 145 
146 - 151 


2 

4 

0 

3 

3 

1 


8*. 15 

16. 17, 18, 19 

10, 12, 14 
11, 13, 20 
9** 


•SUNSPOT CYCLE 8 HAS THE SMALLEST CYCLE DURATION 
EQUAL TO 116 MONTHS 


* ‘SUNSPOT CYCLE 9 HAS THE LARGEST CYCLE DURATION 
EQUAL TO 149 MONTHS 

NOTE : NO CYCLE DURATION HAS BEEN REPORTED WITH A 
VALUE BETWEEN 127 AND 133. 


TABLE 4. STATISTICAL PROPERTIES OF SHORT-PERIOD AND 
LONG-PERIOD SUNSPOT CYCLES 


MEAN 

CYCLE STANDARD 


CYCLE GROUP 

SCN’S DURATION 

DEVIATION 

RANGE 

SHORT-PERIOD 

8. 15-19 

121.8 

3.3 

116-126 

LONG-PERIOD; (1) 

9-14, 20 

140.0 

4.7 

134-149 

12)* 

10-14, 20 

138.5 

3.2 

134-143 

ALL 

8-20 

131.6 

10.0 

116-149 


•EXCLUDES SCN 9 
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A superposed-epoch analysis for each of these two cycle groupings results in 
the curves depicted in Figure 6. It is seen that short-period cycles are characterized 
as cycles v;hich reach maximum more quickly and have a higher than long- 

period cycles. For short-period cycles, R^^ averages about 127 with a standard 


deviation of about 41 (at maximum), reaching maximum in about 44 months. On the 
other hand, long-period cycles average about 92 at R.j, x with a standard deviation 
of about 22, reaching maximum in about 50 months. 



Figure 6. Superposed epoch analysis curves for short-period 

and long-period cycles. 


Returning to Table 1, it is observed that of the six cycles that had cycle 
duration_less than the average cycle duration (i.e., the short -period cycles), four 
had an R^ x value which exceeded the average R MAX value. Also, of the cycles 

with R^ x >. 142, all (3 of 3) have been short-period cycles. Of the seven cycles 

catalogued as long-period cycles, five had an R MAX value lower than average. Thus, 

9 of 13 sunspot cycles fit the pattern where short-period cycles are associated with 
high- valued cycles and long-period cycles with low-valued R^ax c y c ^ es * There- 

fore, a "crude” means is reached whereby cycle duration can be estimated on the 
basis of H max . 

While a plot of R|^ x versus SCN (Figs. 1 and 2) shows a downward trend 

followed by an upward trend, being centered at about SCN 14, it is seen that a plot 
of minimum-to-minimum period versus SCN (Fig. 7) gives the appearance of a "square 
wave," albeit there are really too few sunspot cycles to make a bona fide assessment. 
(See Wilson, et al. [11] concerning the behavior of cycle duration with time.) Place- 
ment of cycle 21 on the plot is, thus, seen to be largely conjectural, being either 
about 140 if the "square-wave" pattern is legitimate, suggesting that cycle 21 is a 
long-period cycle, or being about 122 if the association between short-period cycles 
and high- valued/ cycles is the stronger. If the latter is true, then, either the 

"square-wave" p r " ern is merely illusory or cycle 20 may have been anomalous (see 
Section III.F). 
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Figure 7. Variation of MIN-MIN PERIOD versus SCN. 


E. A Transient Curve-Fit Scheme 

In an earlier paper, Teuber, et al. [2] developed an analytical method for 
fitting impulsive transient phenomena; in particular, for light curves of X-ray 
flares. A formula was derived, given as 

Y = A (Bt) m e“ Bt + YO , (6) 


where Y is the magnitude of the pulse calculated for a particular t, the time from the 
start of the pulse, and YO is an offset or background reference (calculated at t = 0). 
A and B are constants, A being the pulse height and B being the pulse width. 

Values for A and B are computed as 


A = (YM - YO) m" m e m 


(7) 


and 


B = m (ASC) -1 


( 8 ) 


where YM is the maximum magnitude of the pulse, ASC the ascent period of the pulse, 
and m an exponent which had the value of 3/2 in the X-ray flare study. 


Sunspot cycles, perhaps, show an analogous transient behavior and may be 
fitted with equations of the form of equations (6), (7), and (8). Now, YM corres- 
ponds to the maximum smoothed sunspot number at cycle maximum (R^AX^ > YO the 

minimum smoothed 3unspot number at cycle minimum ASC the ascent period of 


the cycle, and t the time in months from cycle minimum. Once these parameters are 
set, m can be varied to determine the best fit. (In practice, one would allow YM, 
YO, ASC, and m to vary to achieve the best fit.) 
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Figure 8 illustrates the effect of varying m. A high value of m means that the 
curve initially begins to change more slowly than lower values of m, but rises more 
steeply as t approaches ASC, Furthermore, a high value of m means that the curve 
falls off more steeply than lower values of m just after peak, but more slowly as t » 
ASC. A fit of the superposed-epoch analysis curves, shown in Figure 6, suggests 
that m has a value of about 3.8 on the ascent side and 4.3 on the descent side for 
short-period cycles, and a value of 3.6 and 3.8 on the ascent and descent sides, 
respectively, for long-period cycles. Thus, using equation (5) to determine an initial 
estimate for and then, assuming that the cycle in question is either a long- 


period or short-period cycle, equations (6), (7), and (8) can be applied to determine 
intermediate values of during the cycle. This initial sunspot cycle curve pre- 
diction can then be adjusted to better fit observed values as they become known. 


(Wilson [5] has found additional ways to deduce R^y 


once the cycle has started. 


These other means are based on the magnitude of the rate of change of smoothed 
sunspot number during its ascent and on the magnitude of the sum of monthly mean 
sunspot numbers, prior to t = 24 months. Both are found to correlate well with 

^MAX^ 



RELATIVE UNIT (NUMBER) 

Figure 8. Variation of transient -fit functional 
shape with value of m. 


F. Application of Empirical Method 

1 . 0 Cycle 21 

Having laid a foundation for estimating sunspot number as a function of time 
from cycle minimum, based on the behavior of cycles 8 through 20, the technique is 
now applied to cycles 21 and 22. The goal, of course, is to more accurately describe 
these cycles and their cycle-related parameters than can be achieved by using mean 
cycle statistics . 
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Based on mean cycle statistics, cycle 21 would be expected to have an 

value of approximately 116,2 ± 36.7 (i.e., the mean value ± 1 unit of standard devia- 
tion accuracy), or 106.9 ± 36.0 if the superposed-epoch analysis curve is used for 
the mean sunspot cycle. It should reach maximum in about 48 ± 6 months, decay to 
subsequent cycle minimum in about 84 ± 10 months, thus having a cycle duration of 
about 132 ± 10 months, and have an of about 53 + 16. If cycle 21 is regarded 

to be a short-period cycle, its R]y| AX should be about 127.1 ± 41.4, its ASC equal to 


about 44 months, and its DES equal to about 78 months, yielding a cycle duration of 
about 122 ± 3 months. Certainly, it would be expected to lie within or very near the 
range of previously observed short-period, cycles which is 116 to 126 months. On the 
other hand, if cycle 21 is a long-period cycle, then its K MAX should be about 91.8 + 


22.0, its ASC equal to about 50 months, and its DES equal to about 90 months, yield- 
ing a cycle duration of about 140 ± 5 months. As a long-period cycle, it would be 
expected to have a cycle duration within or very near the range of previously 
observed long-period cycles which is 134 to 149 months (excluding SCN 9, 134 to 
143 months) . 


Application of equation (5), using the e^x) , 0 2 (x), and 0 3 (x) values contained 


in Table 2, as shown in Figure 2, predicts that cycle 21 should have an Kjyj AX 


of 


about 160 ± 20. This is remarkably close to the value of 154 predicted for cycle 21 
by Sargent [13] and to the actually observed value of 164.5. Using the linear 
regression equations identified in Figures 3, 4, and 5, which are based on compari- 
sons to observed values of Kjyj AX » a SLOPE Agc of about 3.5 ± 0.3, a SLOPE DEg of 

about -1.8 + 0.2, and an E^jeAN of about 71.8 ±5.3 can be projected for cycle 21. 


SLOPE asc for cycle 21 has now been observed and its value is 3.63, very close to 
that predicted from the linear regression equation. Using the observed SLOPE Agc 
value, SLOPEjj Eg for cycle 21 is redetermined to be -1.85. 


ASC for cycle 21 is now known to be 42 months, very close to that expected 
for short-period cycles, and significantly different from the value expected for long- 
period cycles. Likewise, it is now known that for cycle 21 is 164.5, making it 

a high-valued &^ AX cycle, the second highest in the set SCN 8 through 21. It. may 


be surmised, then, that cycle 21 is very probably a short-period cycle. Thus, it 
is expected to have a cycle duration of about 122 ± 3 months, implying that DES = 

80 ± 3 months. This means that for cycle 22 may occur as early as August 

1986 ± 3 months, much earlier than the June 1987 date calculated from the average 
11-year cycle duration, or later dates if SCN 21 is a long-period cycle. Yoshimura [14] 
has also suggested that cycle 21 will be a short-period cycle, as have Otaola and 
Zenteno [15]. Yoshimura bases his prediction on a repeating pattern of "grand 
cycles," each of 5-cycle duration; Otaola and Zenteno base their prediction on appli- 
cation of a non-integer technique of power spectral analysis. Note that if cycle 21 
is indeed a short-period cycle, then its minimum-to-minimum period will not fit the 
supposed "square-wave" pattern depicted in Figure 7, implying, perhaps, that the 
pattern is merely coincidental or that cycle 20 may have been anomalous. When cycle 
20 is compared to the most recent cycles, it certainly appears to have been anoma- 
lous, especially in terms of its Rjyj AX value. (See Wilson, et al. [11] for additional 

comments concerning the behavior of cycle duration with time.) In this cited paper, 
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based on an investigation of "trends" in the sunspot record, a statistic is identified 
which strongly suggests that cycle 21 is a long-period cycle, with for cycle 22 

occurring about February 1988, very close to that predicted by Holland and Vaughan 
[16], based on a Lagrangian least-squares fit of smoothed 2800 MHz radio flux. (See 
also Xu, et al. [17] .) 


In addition, it must be noted that Fairbridge and Hameed [18] have predicted 
cycle 22 minimum occurrence to be 1989.1 ± 0.9, on the basis of a suggested repeat- 
ing pattern of sunspot minima every 178 years; thereby, implying that cycle 21 is a 
long-period cycle of period equal to about 151 months, almost 4 standard deviation 
units longer than the mean for long-period cycles (excluding SCN 9). (See also 
Jose [19] and Gregg [20].) A period of 151 months implies that the descent period 
for cycle 21 will be 109 months, nearly 2.6 standard deviation units longer than the 
mean descent period for cycles 8 through 20. Together, this seems to suggest that 
cycle 21 is very improbably a long-period cycle, being much more likely a short- 
period cycle. The analysis of Fairbridge and Hameed, as that of Jose, Gregg and 
the other aforementioned analyses, is based, in part, on sunspot data prior to SCN 8. 
It is recalled that sunspot data prior to SCN 8 are not considered to be as reliable as 
that for cycles occurring after SCN 8. This may mean that the supposed association 
between sunspot minimum occurrences and the 178-year rhythm of the solar orbit 
around the center of mass of the solar system is not as strong as has been previously 
believed. (Note, however, that if the lower limit of the Fairbridge and Hameed 
result is used, namely that SCN 22 will begin about 1988.2, then this is consistent 
with SCN 21 being a long-period cycle with cycle minimum for SCN 22 occurring 
about February 1988.) 


Figure 9 compares observed Et 13 values with those predicted for cycle 21, based 
on the prediction of equation (5) and the observed S MIN value. Two SLOPE Agc 

and SLOPEjjgg lines are drawn. The heavy lines are associated with the observed 
SLOPE Agc value; the light lines are associated with the predicted SLOPE Agc value. 



o 

JUN 1*7« 


20 » SO 100 120 

I TIME IN MONTHS MOM W MIN OCCURRENCE 


Figure 9. Comparison of observed and predicted smoothed sunspot numbers 
for SCN 21 using SLOPE Ag( -, and SLOPE^gg. 
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The dots represent observed values. The figure shows that SLOPE^qq, whether 

the predicted slope or the observed slope is used, fits the observed values 

reasonably well, especially for t > 20. Cycle 21 is observed to have had a period of 
about 6 months, just after R MIN occurrence, when K lg remained essentially constant. 

At t = 6, observed Rj, g values are seen to rise nearly linearly to occurrence, 


which was at t = 42. There has been one significant "bump" in the decline of cycle 
21, reminiscent of the behavior of cycle 20 (which was a long-period cycle), occurring 
approximately 21 months after S MAX occurrence. 


Figure 10 compares a fit, based on equations (6), (7), and (8), with the 
observed R^ g values for cycle 21. The comparison is made for two different epochs: 

(1) post-H^j^ occurrence, but prior to Rj^ ax occurrence and (2) post 

Sjviax occurrences. For the first case (Fig. 10, bottom), = 12.2 and * s 

assumed to be 160, based on the prediction of equation (5). An m equal to 3.8 on 
the ascent side and 4.3 on the descent side is used; i.e., the mean m values for the 



CASE l! Bmin - 12,2, %AX " 160 ' Asc ' 44 - 
200 r m ASC * 3.8, m 0ES • 4,3 



t. TIME FROM R M ( N OCCURRENCE 

Figure 10. Comparison of observed and predicted smoothed sunspot numbers 
for SCN 21 using the transient-fit scheme for two selected cases. 
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short-period cycle group have been applied. An ASC of 44 months is employed. It 
is noted that there is very close agreement between observed and predicted values 
for 0 < t £ 55. For t > 55, corresponding to the "bump" in the decline of cycle 21, 
the curves are in disagreement. Observed values of K 13 have remained significantly 

higher than the plotted curve. It is as though the predicted decline curve has been 
shifted about six months to the right. The "bump" can also be approximated os the 
superposition of two "transient" pulses, with peaks separated by about 1.5 to 2 
years, although such a fit is not shown here. This may be an indication that this 
cycle, like its predecessor cycle 20, may have one or more "bumps" during its 
decline and that these "bumps" may be the result of short-term periodicities in the 
sunspot cycle, shorter than the 11-year average cycle duration (see Wilson [5,6] and 
references contained therein). 

For the second case (Fig. 10, top), = 12.2, = 164.5, ASC = 42, 

and m = 2.9 (the "best-fit value of m for 0 ^ t < 78; i.e., data through December 
1982). The curve shown represents the "best-fit" using all data from June 1976 
through December 1982, which includes the "bump." A "best-fit" for all data prior 
to the "bump," i.e., 0 ^ t ^ 55, yields a value for m of about 5.0. Thus, the effect 
of the "bump" is to significantly lower the value of m during descent, at least for 
t < 78. If cycle 21 is a short-period cycle, the "best-fit" m value must increase in 
value for t > 78, which effectively lowers the predicted R l3 values. Otherwise, at 

t = 122, a time when R 13 should be at or very near the predicted R 13 values 

will be too high. For m - 5.0, the standard deviation s about the predicted curve 
(not shown) for the ascent period t s 42, is about 5.4. For m = 2.9, the "best-fit" 
for t ^ 78 including the "bump," s is about 13.0. 


2.0 Cycle 22 


Cycle 22 begins when R^ 3 has reached a minimum value in the decline of cycle 

21. As stated above, based on the statistics of cycles 8 through 20, cycle 21 is 
expected to be a short-period cycle; consequently, its minim urn -to-minimum period 
should be about 122 ± 3 months. As a short-period cycle, cycle 21 Would be expected 
to be less than 128 months and more than 116 months in duration. Hence, R^j^ for 

cycle 22 should occur between February 1986 and February 1987, with the "best- 
guess" being about_August 1986. If cycle 21 is a long-period cycle, as suggested 
by Figure 7, then R^^ for cycle 22 will occur about February 1988. 


Application of equation (5), using the 6 1 (x), 0 2 (x), and 0 3 (x) values for 


cycle 22 in Table 2, predicts that cycle 22 will have an 


Ejyjax Vfl i ue °f about 110 ± 


20 (Fig. 2), a value being about that observed for SON 20. Using the linear regres- 
sion equations identified in Figures 3, 4, and 5, which are based on comparisons to 
H Max » we Project SLOPE Agc to be about 2.2 ± 0.3, SLOPE DEg -1.3 ± 0.2, and 

R MEAN ^* 8 ^ 5,8t A value of Rj^ ax = i m PR es that cycle 22 will be a low-valued 

R max cycle; consequently, it is expected to be a long-period cycle. As a long- 


period cycle, it should have an ASC of about 50 months and a total cycle duration of 
about 140 months. 
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On the basis that cycle 22 will be a long-period cycle with ASC = 50 months, 
Figure 11 (bottom) plots the anticipated SLOPE^g^ and SLOPE^^g lines. This 

approach yielded quite reasonable values for cycle 21 (recall Fig. 9), so it is antici- 
pated that it will yield reasonable values for here. Therefore, cycle 22 is pre- 
dicted to reach a maximum of = HO ± 20 about October 1990 (±6 months), and 

to reach subsequent cycle about April 1998. (These dates are based on cycle 

21 being a short-period cycle. Other dates have been summarized in the abstract, 
on the basis of SON 21 and 22 being either long-period or short-period cycles.) A 
projection for cycle 22 using the "transient 1 ' fit approach is also depicted in Figure 11 
(top), where = 5, Et MAX = 110, ASC = 50, and m - 3.7, the value for m midway 

between m AS( -, and m^g for long-period cycles. It is emphasized that Figure 11 

presents the "best guess" for the general behavior of cycle 22. Until cycle 22 com- 
mencement has occurred, this prediction cannot be improved. Also, the technique 
in its present form does not allow the anticipation of the occurrence of "bumps" 
following occurrence, as has been seen for cycles 20 and 21. 


transient fit 

Emin * &. Wmax * i io. asc • m m » 3 7 




Figure 11. Prediction for SCN 22 using the SLOPE. and SLOPE ni7Q 

fit and the transient fit. 
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Another estimation for cycle 22, on the basis of paired-cycle duration using 1 
principal component analysis, is that cycle 22 H«, AX will be about 120 ± 30 (Teuber 
[12]). This is consistent with what has been predicted here, and also to that pre- 
dicted by Holland and Vaughan [16] based on a Lagrangian least-squares prediction 
of smoothed 2800-MHz radio flux. They predict cycle 22 to have a peak smoothed 
2800-MHz radio flux of about 150 Which, based on a linear regression comparison 
given by Wilson [5], yields a peak smoothed sunspot number of about 95. Also, 
Xu, et al. [17] have predicted cycle 22 to have an between 90 and 106, It 


must be noted, however, that these predictions are in disagreement with that pre- 
dicted by Gregg [20], based on a novel oscillator model; he predicts cycle 22 to 
begin in 1988, peak about 1994 with an annual mean sunspot number of about 36, and 
have a descent period shorter than its ascent period, reaching the following minimum 
in 1999. (Wilson [6] has shown that maximum annual values of sunspot number are 
very close in value to K MAX values; for a maximum annual value of 36, & MAX is 


expected to be about 39.) Negative skewness, that is ascent duration being longer 
than descent duration, is non-existent in the sunspot data for cycles 8 through 21. 

If one includes cycles 1 through 7, as has been done by Gregg, then there may have 
been at least a few examples of negative skewness. (We remind the reader that 
sunspot data for cycles prior to SCN 8 are considered to be unreliable.) 


IV. CONCLUSIONS 


In this study, a simple empirical method has been developed to predict certain 
parameters of future solar activity cycles — e.g., R max , ASC, DES, MIN-MIN 

PERIOD, SLOPE asc , and SLOPE DEg . The method is based on cycles 8 

through 20 and incorporates a 3-component sinusoidal fit of Rj^ ax values versus 

sunspot cycle number, where the three sinusoidal components include a 90-cycle, an 
11-cycle, and a 2-cycle periodicity. (Linear regression equations, based on 

allow for the evaluation of the aforementioned cycle-related parameters.) A "transient* 1 
curve-fit scheme is also presented which allows for a month-by-month prediction of 
smoothed sunspot number from occurrence for a given cycle. The empirical 


method has been applied to cycle 21, as an example, and found to work very well. 

value of 


Had it been used early in the cycle, it would have predicted an R 
about 160 for cycle 21, compared to the actually observed Hjyj AX 


MAX 

value of 164.5. i If 


cycle 21 turns out to be a cycle of short duration (that is, a cycle having a duration 
of 122 ± 3 months) rather than a cycle of long duration (duration 140 ± 5 months), 
then cycle 22 should begin about August 1986, certainly before February 1987. An 
5*MAX value of 110 ± 20 is projected for cycle 22, occurring about 50 months after 

Rj^IN occurrence, or about October 1990 (+6 months). Cycle 22 is projected to be a 

cycle of long duration; hence, it should end about April 1998. If instead cycle 21 
turns out to be a cycle of long duration, then cycle 22 will not begin until about 
February 1988, peak until about April 1992 and end until about October 1999, 
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